Geometry of depolarizing channels 
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Depolarizing maps acting on an TV dimensional system are completely positive maps resulting into 
compression of the Bloch 'ball' along N 2 — 1 polarization directions. In the qubit case these maps 
are a convex sum of four extremal maps and form a simplex in the space of compression coefficients 
along the three polarization directions. We calculate the compression domain for three and four 
level systems. For a three level system the region has curved surfaces, but it is a simplex for a four 
level system. We conjecture that it is a simplex in the case of 2" level systems. 

PACS numbers: 



I. INTRODUCTION 

The dynamics of open quantum systems can be de- 
scribed in terms of Completely PositivefCP), trace pre- 
serving maps acting on the system [1, U, B S S @| ■ In 
general these maps can rotate, compress and translate 
the Bloch ball, resulting into 'Affine' maps [HI]. These 
maps are used to simulate different noises that could be 
acting on the system [l|, 0] ■ 

A simple example of noise can be given by depolarizing 
maps, which result into shrinking of the Bloch ball. In 
general the compression can be anisotropic, giving a rich 
geometric structure to the problem. The requirement of 
complete positivity turns out to be stronger than just 
positivity. For example, for N=2 the map corresponding 
to shrinking of the Bloch sphere only along one direction 
and leaving the other two invariant (pancake map) is not 
CP. In N=2 case the general depolarizing map has a very 
simple structure and can be written as a convex sum of 
four fixed extremal maps. One would like to see if it true 
for higher dimensional systems. 

In what follows we answer the question for N=3 and 
N=4 cases. We also give a way to answer the question 
for higher level systems, although we must admit that 
the calculations becomes much more involved. 

Our paper is organized as follows: In Sec. II we briefly 
describe the geometry of N=2 level depolarizing chan- 
nels. In Sec. Ill we show how to calculate the region of 
positivity for an N dimensional system. We carry out the 
calculations for N=3 and N=4 cases in Sec. IV and V 
respectively. We conclude the paper with further discus- 
sions in Sec. VI. 



II. THE QUBIT CASE 

The qubit density matrix is given by 



Positivity of p requires that a lies inside a sphere with 
unit radius. Let's assume that due to depolarization, the 
sphere shrinks along the three polarization directions by 
factors of v\, v 2 and v 3 . The natural question arises - 
what are the allowed values of v\, v 2 and v 3 to ensure 
complete positivity of the map? 

It turns out [l2j that the allowed region in the space of 
v\, V2 and ^3 forms a tetrahedron with vertices at (1,1,1), 
(1,-1,-1), (-1,1,-1) and (-1,-1,1). These vertices corre- 
spond to four extremal maps given by Bi(p) — utpai, 
where i = 0, 1, 2, 3 and uq = I 2 - The general depolariz- 
ing map can be written as a convex sum of these maps 
as 



B(p) = paipai + qa 2 po- 2 + ra 3 pa 3 + (l-p-q-r)p (2) 



where p, q, r are non-negative. 

Only the region inside the tetrahedron is CP. For ex- 
ample, partial transpose corresponds to the point (1,-1,1) 
and lies outside the region. 



p = ~ (7 + a ■ a) 



(1) 




FIG. 1: As a result of depolarization the Bloch sphere shrinks 
anisotropically 
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FIG. 2: CP region in the space of i/'s is a tetrahedron with 
vertices corresponding to maps lpl, <j\po\, oipo-z and ozpoz 

III. THE N LEVEL SYSTEM 

For an N level system the density matrix can be writ- 
ten in terms of N 2 — 1 traceless, Hermitian and trace- 
orthogonal generators of SU(N), denoted by Ji, and the 
polarizations along those directions (aj) as 

This normalization results into pure states having unit 
length [JJJ], analogous to the qubit case. 

Let us assume that a general map acting on p takes it 
into another density matrix given by p' . The map can be 
written in terms of the A matrix [3| as: 



For completeness we take T to be TV 2 g) N 2 dimensional 
matrix acting on the vector (1,7?). Needless to say, T\\ — 
1 as required by the trace preserving condition. Both 
the T and the A matrices contain all the information 
about the map and are thus related. The relation can be 
established in the following way. 

Since the super-matrix A acts on all the entries in p, 
we can convert p into a vector (let's call it V) which is 
acted upon by A. The entries in V are linearly related 
to the coherence vector as 

VP = M ijaj (8) 

The map takes V p to V which is related to V p as 

V p ' = M ■ a' = M ■ T ■ a = A ■ V = A ■ M ■ a. (9) 

Thus 

A = M-T-M- 1 . (10) 

Once we get A and B we can find the eigenvalues of B 
and check for complete positivity. We solve the case of 
N=3 and N=4 in next two sections. 



IV. N=3 CASE 

The qutrit density matrix is written as 

1 8 
p= (l + V3^aA) (11) 
6 i=i 

where X's are the Gell-Mann matrices 



Prs A rsr i s ip r 



(4) 



With change of indices, we can find another matrix B 
which is Hermitian and has simpler conditions of com- 
plete positivity. It is obtained by 

B rr ' ss ' — A rs r ' s i. (5) 

Complete positivity requires B to be positive semi- 
definite (and thus have non- negative eigenvalues). 



A. Deriving A from compression of the coherence 
vector 

The simplest way to understand the depolarizing map 
is to think of it as compression of the coherence vector. 
Let's assume that as the result of the map the polariza- 
tions are modified as 



a l = Vidi. 

which can be considered as a special case of 



(6) 
(7) 



Ai 



At 




As shown in ([7]), T can be written as 
T = diag(l, v\,V2, —, Ms) 



1 
-1 



1 



(12) 



(13) 



M relates the density matrix to the coherence vector 



as: 



VP = .1 /,//,. 



(14) 



Explicitly: 
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Calculation for A = M ■ T ■ M" 1 yields: 
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The B matrix is obtained by change ol indices as 



and is given by: 



(17) 



B 



/ 2+3v 3 +VS 

' 6 






V± I V2 

2^2 





\ 2 "I" 2 





2 — 3f3+"8 



ill 

2 



£2 

2 











vi 1 v 2 
2 ' 2 























V4.+V5 

2 


\ 







2 2 

































vg 
3 















2 




V5 

2 



















2— 31/3 + ^8 

6 







































2+3v 3 +V8 

6 























V&+V7 

2 















1 
3 




V8 

3 


1 

3 







v e 
2 




v 7 
2 









V5 

2 

















vs 
3 


1 

3 























v§ 
2 




"7 

2 











V8 

3 















2^2 























l+2v s 
3 


/ 



(18) 



For the map to be positive, all the eigenvalues should 
be positive semi-definite. Six of these eigenvalues are 
given by the 'hyperplanes': 
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The remaining three eigenvalues are given by the eigen- 
values of 



2+3^3+^8 
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(20) 



Since the matrix is Hermitian and thus has real roots, 



the condition of positive eigenvalues translates to (A7 + 
Ag + Ag),(A7Ag + AgAg + AgA7),A7A s Ag being positive. 
These can be easily found from the characteristic equa- 
tion as 



h-j — A7 + Ag + Ag 

= 1 + ^3 + i/ 8 (21) 

■Sl = A7A8 + AgAg + AgA7 

= ^ (~ 3 ^ + V2) 2 + 3 ^ - 3 + ^5) 2 - 3 (ye + vt) 2 + (2 + 1*) (2 + 3i*) + 2u z (4 + 5i*)) . (22) 

S2 = A7AgAg 

= -j-(54^ 4 + 1/5) K + vi) - 18^(1 + 2^ 8 ) - 18i/|(l + 2v & ) + 18i^(3(i/ 4 + ^)(>6 + *y) - 2i*(l + 2i/ 8 )) 

+ (2 + 3^3 + i*)(4 - 9(^4 + ^ 5 ) 2 - 9(i/ 6 + ^ 7 ) 2 + 10^ 8 + 4^ 8 2 + 6^3(1 + 2u s ))). (23) 

I 



Hence we see that the completely positive region is 
bounded by seven hyperplanes, a second degree surface 
si and a third degree surface S2. This is quite different 
from the qubit case where the allowed region is bounded 
by four planes and is just a simplex with four vertices. 
For N=3 the depolarizing map has an infinite number of 
extremal points corresponding to the curved surface. 

V. N=4 CASE 

From the previous section we might expect that the 
CP region will be a curved surface for N=4, but quite 



analogous to the qubit case it turns out to be a simplex. 

Let's start with taking the generators to be: J={I®ai, 
(Xi®I, a\®CTi, cj2®(Ji, a?,®(Ji \ for i=l, 2, 3. The density 
matrix is written as [3, [lCj, LLl| 



p = - (I + Vda ■ J) (24) 



The M matrix is given by: 
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Let's take the T matrix as in the previous case B matrix is given by 

T = diag{l,v 1 ...v 1B }. (26) 
I 
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The eigenvalues of B are 
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Since all of the eigenvalues are hyperplanes, CP region 
is a simplex. Also, the vertices of this polyhedron cor- 
respond to the maps lpl, JipJi, identical to the qubit 
maps. The extremal map - JtpJi is just a rotation of 
the Bloch ball about the J, axis. The general depolariz- 
ing map can be written as a convex sum of these sixteen 
extremal maps. 

VI. DISCUSSION 

We found that the general depolarizing map in N=2 
and N=4 cases has a simple geometric structure, while 
for N=3 it is much more complicated. We conjecture 
that for any 2 ra dimensional case the CP region will be a 
simplex, with vertices corresponding to a unitary rotation 



fa 


+ ^9 


- v w 


- V\\ 


- V\i 


- Vl3 


- ^14 


- v 15 ) 


V8 


- v% 


+ v w 


+ V\\ 


+ Vyi 


- V13 


- Vx4 


- v 15 ) 


^8 


- v% 


+ vio 


- v\\ 


- V\1 


+ "13 


~ Vx4 


- V15) 


^8 


+ V<A 


- v w 


+ V\\ 


+ Vyi 


+ Vl3 


- VX4 


~ V\h) 


v$ 


- v% 


- v w 


+ V\\ 


- V\1 


- "13 


+ V\4 


- ^15) 


^8 


+ ^9 


+ v w 


- V\\ 


+ Vyi 


- V13 


+ via 


- V15) 




+ Vg 


+ v w 


+ V\\ 


- V\1 


+ Vl3 


+ VxA 


- v 15 ) 


v$ 


- v% 


- v w 


- V\\ 


+ Vxi 


+ V13 


+ via 


- V15) 


^8 


- v% 


+ v w 


+ V\\ 


- V\i 


- "13 


- Vli 


+ vi 5 ) 


^8 


+ Vg 


- v w 


- V\\ 


+ VX2 


- "13 


- VxA 


+ V15) 


^8 


+ Vg 


- v w 


+ V\\ 


- VX2 


+ Vl3 


- VxA 


+ vi 5 ) 


VB 


- Vg 


+ v w 


- V\\ 


+ VX2 


+ V13 


- VxA 


+ vi 5 ) 


v$ 


+ Vg 


+ v w 


- V\\ 


- VX2 


- Vl3 


+ VxA 


+ vi 5 ) 


v$ 


- Vg 


- v w 


+ V\\ 


+ Vx2 


- "13 


+ via 


+ vi 5 ) 


^8 


- Vg 


- vio 


- v\\ 


- VX2 


+ "13 


+ VX4 


+ v 15 ) 


^8 


+ ^9 


+ vio 


+- V\l 


+ V\2 


+ fl3 


+ VXA 


+ v 15 )) 



(28) 
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of the Block ball. It would be interesting to explore the 
extent to which the N=2 and N=4 cases are similar for 
other maps. 

We would also like to point out that our method of 
calculating B, and thus finding out the region where a 
map is CP is general and is not restricted to depolarizing 
maps. The most general map (Affine map) can be easily 
obtained if we let T also contain translations (see [12|). 
In that case T becomes 

where T' contains only the unital part. 
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